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Abstract 

A condition for the synchronizability of a pair of PDE systems, coupled through a finite set of 
variables, is commonly the existence of internal synchronization or internal coherence in each sys- 
tem separately. The condition was previously illustrated in a forced-dissipative system, and is here 
extended to Hamiltonian systems, using an example from particle physics. Full synchronization is 
precluded by Liouville's theorem. A form of synchronization weaker than "measure synchroniza- 
tion" is manifest as the positional coincidence of coherent oscillations ("breathers" or "oscillons") 
in a pair of coupled scalar field models in an expanding universe with a nonlinear potential, and 
does not occur with a variant of the model that does not exhibit oscillons. 

PACS numbers: 05.45.Xt,05.45.Jn,11.27.+d 

Keywords: chaos synchronization, coherent structures, measure synchronization, oscillons, scale interactions 



1 (May 30, 2009) 



The phenomenon of synchronized chaos, initially explored in low-order ODE systems 
, has been extended to PDEs that describe a variety of systems of physical interest 
Chaos synchronization_extends the paradigm of synchronization of regular oscillators that is 
ubiquitous in Nature [5] . One seeks an understanding of the internal properties of a chaotic 
physical system that will allow a pair of such systems, loosely coupled, to synchronize, 
despite sensitive dependence on initial conditions. Spatially extended systems offer richer 
possibilities for relationships that fall short of full synchronization than do ODE systems. In 
geophysical example previously studied [gI, it was seen that slaving of small scales, a relation- 
ship that defines an inertial manifold, was crucial to synchronizability. Here, it is suggested 
that more general inter-scale relationships, as may give rise to coherent structures within 
each system separately, are required for the synchronizability of the pair. The connection is 
illustrated in a particle physics context - a toy model of a scalar field in the expanding early 
universe a Hamiltonian system without an attractor. Weak scale relationships allow 
oscillons (breathers) to persist, and give rise to a new form of synchronization defined by 
the coincidence in position of oscillons in two coupled scalar field models. Conversely, where 
the dynamics do not admit such coherent structures, vestiges of synchronization are lost. 

In the non-Hamiltonian meteorological example described in which we first review, 
two geophysical fiuid systems, representing planetary- scale wind patterns, are coupled only 
through their medium-scale Fourier components. Each system is given by a potential vor- 
ticity equation 

^^^ + Ji^P.,q^) = F, + D, (1) 

where the streamfunction ip is the fundamental dynamical variable, the Jacobian 
J(^/', ■)=^^ — gives the advective contribution to the Lagrangian derivative D/Dt, 
there are two horizontal layers i=l, 2, and the potential vorticity o, which generalizes an- 
gular momentum, is a derived variable defined in terms of in 61]. Potential vorticity is 
conserved on a moving parcel, except for forcing Fi and dissipation Di. 

Two models of the form ([I]), Dq^/Dt=F^+D^ and Dq^ /Dt=F^+D^ were coupled 
diffusively through one of the forcing terms: 



Pi = + (2) 

— * 

where the flow has been decomposed spectrally and the subscript k on each quantity indicates 
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the wave number k spectral component (suppressing the index i). A background flow q* 
forces each system separately. The set of coefficients /i| was chosen to couple the two 
channels only in some medium range of wavenumbers. Band-limited coupling defined by /i| 
replaces the coupling of two PDE systems at a discrete set of points as in [4|. 

The two systems synchronize over time in Fig. [H where the contours of ip are streamlines 
that define the flow. The synchronization is manifest as the coincidence of structures of 
meteorological significance ("blocking patterns" that interrupt the flow) in both space and 
time. That coincidence is robust against significant differences in the two systems jsl. 

The large scales need not be coupled because of well known dynamical relations between 
scales ( "inverse cascade" ) in 2D turbulence 9| . But the stronger result 6| that synchroniza- 
tion occurs without coupling of the smallest scales is explained simply: The smallest- scale 
components are thought not to be independent dynamical variables but are functions of 
the medium and large-scale components, defining a dynamically invariant inertial manifold. 
Where the set of variables that are coupled is smaller than a minimal independent set, 
synchronization does not occur. Partial synchronization of forced-dissipative systems will 
usually be easy to arrange, as previously found by Kocarev et al. jj], since approximate 
inertial manifolds (AIM's) exist for almost all parabolic PDEs lOj. 

In Hamiltonian systems, the subject of the current work, inertial manifolds cannot exist, 
since the stability of such a manifold would imply a collapse of phase space volumes for 
trajectories that start off the manifold, contradicting Liouville's Theorem. A simple example 
is the Klein- Gordon equation in an expanding background geometry, in one space and one 
time dimension, with cyclic boundary conditions, with a nonlinear potential of a type that 
gives rise to oscillons, possibly representing the first coherent structures in the universe. The 
field satisfies 

where if is a Hubble constant, the potential V is given by V{(j)) = (1/2)0^ — (1/4)0"^ + 
(1/6)0^ for units in which the scalar field mass m = 1, and the prime denotes the derivative 
with respect to 0. The equation ([3]) is derived by replacing derivatives in the standard 
Klein-Gordon equation by covariant derivatives for a Robert son- Walker metric describing 
expansion with Hubble constant H llj, giving a term in dip/dt which formally resembles 
friction but which in fact preserves the Hamiltonian structure. Indeed, the system ([3]) is 
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derivable from a time-dependent Hamiltonian density 

n = (l/2)e-^*[(0.)2 + Tc'] + e^V(0) (4) 

where vr = e^*0 is the canonical momentum that is conjugate to 0. Liouville's theorem 
applies even with the time dependence, so that neither an inertial manifold, in the strict 
sense, nor an AIM with a usefully small approximation bound can exist. 

Oscillating coherent structures appear as regions of high energy density at final time, in 
the numerical integration shown in Fig. [2^, after initializing with thermal noise as in [3]. 

Consider two systems ([3]) coupled diffusively, through some Fourier components of the 
field only, according to 

ff±A,B rlfh"^'^ rPrh"^'^ 

with the coupling coefficients Ck = vanishing above a threshold value \k\ > ko and set to 
a large value for \k\ < ko, so that corresponding large scale components in the two systems 
are effectively clamped. For coupling of the lowest wavenumber modes, up to wavelengths 
of about twice the final oscillon width, oscillons in the two systems were found to occur 
at mostly the same locations, though their amplitudes differed, as shown in Fig. [3^. In 
contrast, if the two systems were left completely uncoupled, but shared initial conditions over 
a range of scales, and only the smallest scales, well below the oscillon widths, were initialized 
differently, then oscillons formed at locations that were uncorrelated (Fig. [3]d). Apparently, 
there is a "butterfly effect" , as in meteorology, through which the small scales have a large 



impact on the positions of formation of coherent structures. [12] But the dynamical evolution 
of these structures then partially slaves a portion of the small-scale sector (as with shock 
waves) and proceeds independently of the remaining portion. 

With severely attenuated initial noise (Fig. H^), oscillon synchronization occurs with an 
even narrower range of coupled wavenumbers. The large-scale components (top 32 modes) of 
the same oscillons are shown in the bottom portion of Fig. |Dd. For four of the five oscillons 
shown, the oscillon positions are coincident with the local maxima of the truncated field. 
For the remaining one (the leftmost), the oscillon is displaced from the local maximum 
of the truncated field, but is still coincident with the corresponding oscillon in the other 
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system (FigH^). The small-scale Fourier components, in both situations, are slaved to the 
large-scale components insofar as they determine oscillon position. 

Coincidence of oscillon positions suggests measure synchronization, the weak form of 
synchronized chaos in which the trajectories of two coupled systems become the same when 
the systems are coupled, without a requirement that the states of the systems are the 
same at a given instant of timejisl. [l^. Measure synchronization is characteristic of jointly 
Hamiltonian coupled systems. Here, the coupled oscillon systems do not quite attain measure 
synchronization, since the corresponding oscillons differ in amplitude, and the configuration 
([5]) which can be written: 

^f'^ = e-"W = dn^'''/d4''' (7) 

■ A.B ^n-/A,BipiiA.B . / B,A A,B\ /o\ 

is not derivable from a joint Hamiltonian. (It would be Hamiltonian if the second term in ([8]) 
were to be replaced by Ck{4>^'^ — (pt'^)-) -^^^ even without the joint Hamiltonian structure, 
complete synchronization is not likely achievable by coupling a finite number of modes, if 
the remaining infinite number are not slaved. 

In a comparable system without oscillons, there are no vestiges of synchronization. Cor- 
relations between corresponding modes in the coupled systems are displayed in Table [TH 
for the oscillon system, and in Table [T)d for an alternate pair of systems with potential 
V = Vg^]^^{(j)) = (1/2)0^ + (1/4)0'' -|- (1/6)0^ in ([3]), for which no oscillons form, as seen 
in Fig. [2]d. (Oscillons occur in one case and not the other because the —0^ term in the 
first case gives a flatter potential, so that larger amplitude oscillations have lower frequen- 
cies and decouple from the faster, smaller travelling wave solutions that would otherwise 
cause them to dissipate.) While there are small but sigificant correlations between some 
of the uncoupled modes in Table [IH, corresponding to the small partially-slaved portion of 
the small-scale sector, no significant correlations appear between corresponding small-scale 
components that are not coupled in Table [Do. The correlations in the case with oscillons are 
expected to extend to much shorter wavelengths in longer simulations that approach the os- 
cillon lifetimes, as the oscillons decrease in width (in comoving coordinates) and background 
fluctuations decrease in amplitude 

The synchronization results for oscillons suggest a unifying principle governing synchro- 
nization in both Hamiltonian and forced-dissipative systems: Two complex systems of either 
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type can be made to synchronize, with a restricted set of coupled variables, if and only if 
each system exhibits synchronization internally. Sufficiency of the condition follows from 
transitivity if the internal synchronization is exact: Any external coupling that does not 
destroy the internal relationships and causes some pairs of variables in the two systems to 
agree will also cause all variables that internally synchronize with the coupled variables to 
agree with their counterparts in the other system. In the cases considered here, the condi- 
tion also appears necessary. Localized coherent structures like oscillons, as well as long-range 
connections such as the Atlantic-Pacific "teleconnection" described previously , can be 
regarded as approximate forms of internal synchronization. Where a master set of variables, 
in Fourier space, synchronizes with the entire remaining infinite set, an inertial manifold 
exists. Local coherence has been described as synchronization of a more limited type in 
PDEs [3] and in lattices of coupled maps 



The internal synchronization condition also provides guidance as to the choice of coupled 
variables. If the extended systems are "truth" and a "computer model" to which truth is 

n 

coupled in one direction only, via observations, as in meteorology [17[ , then a set of variables 
should be observed that is synchronized with a maximal set of other variables internally. 

An internal coherence criterion applies to ODE systems, consistently with Pecora and 
Carroll's criterion of negative "conditional Lyapunov exponents" in the Lorenz system 
Lorenz X can be said to approximately synchronize with Lorenz Y, along the near-planar 
Lorenz attractor, in agreement with the well known fact that sufficiently strong coupling of 
either X or Y in one system to the corresponding variable of a second Lorenz system will 
synchronize the two systems. Coupling of the Z variables, that do not correlate internally 
with either X or Y, will not do so. But the internal coherence criterion is much more 
useful for PDEs that describe spatially extended systems, both because such systems are 
less tractable analytically - a full set of conditional Lyapunov exponents is hard to compute 
- and because the coherent structures are more meaningful physically. 

For Hamiltonian systems coherent structures seem to bear strongly on synchronizability 
and on the form of synchronization. The ubiquity of coherent structures in solutions to 
nonlinear PDE's is the basis on which potentially broad relevance is claimed for the current 
work. Solitons, for instance, in a pair of PDF systems that exhibit them, might also be made 
to coincide in position, and to move in synchrony, if only a restricted set of corresponding 
Fourier components of the two fields are coupled. The question is whether the internal 
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dynamics that allow the structures to exist in each system would cause the needed parts of 
the uncoupled components to follow. The phenomenon, if it exists, might be used for secure 
communications in the same manner as the systems considered in ^. 

The present generalization differs from measure synchronization [isl in several notable 
respects: First, in an ergodic system with trajectories that define a uniform measure, such 
as the system with modified potential V^^^t ^^^^ exhibits no oscillon behavior, measure syn- 
chronization is trivial. Second, as already pointed out, positional coincidence of coherent 
structures that differ in amplitude or detailed shape is more general than measure synchro- 
nization. Third, there is no requirement that the combined system be Hamiltonian to exhibit 
the weak form of synchronization described here. The impossibility of strict synchronization 
follows not from the joint Hamiltonian structure, but from the impossibility of slaving or 
approximately slaving all uncoupled variables in each system separately. 

The new types of synchronization appear to be equivalent, in their details, to detailed 
inter-scale relationships, within each system internally, that in the present case allow oscillons 
to persist and to stably maintain their positions. The partial agreement of small scales 
between the two synchronized systems, sufficient to force coincidence in oscillon position, 
when large scales are clamped, indeed defines a partial slaving that counters the butterfly 
effect. For consistency with Liouville's Theorem, there must be a compensating expansion 
of the remaining part of the uncoupled modes in phase space, as entropy is cast off to scales 
that are yet smaller. Details of the partial slaving remain to be worked out. The existence 
and the form of synchronization may provide diagnostics for the slaving relationships. 

Acknowledgements: This work was supported in part by NSF Grant 0327929. The author 
thanks Joe Tribbia, Roger Temam, Noah Graham, and Jeff Weiss for useful discussions, and 
especially thanks Alan Guth for also conducting some of the numerical experiments. 



[1] L.M. Pecora and T.L. Carroll, Phys. Rev. Lett. 64, 821 (1990). 

[2] V.S. Afraimovich, N.N. Verichev, and M.I. Rabinovich, Inv. VUZ Radiofiz. RPQAEC 29, 795 
(1986). 

[3] H. Pujisaka and T. Yamada, Prog. Theor. Phys. 69, 32 (1983). 

[4] L. Kocarev, Z. Tasev, and U. Parlitz, Phys. Rev. Lett. 79, 51 (1997). 



7 (May 30, 2009) 



[5] S. Strogatz, Sync: The Emerging Science of Spontaneous Order, Theia, New York (2003). 
[6] G.S. Duane and J.J. Tribbia, J. Atmos. Sci. 61, 2149 (2004). 

[7] E. Farhi, N. Graham, A.H. Guth, N. Iqbal, R.R. Resales, and N. Stamatopoulos, Phys. Rev. 

D 77, 085019 (2008). 
[8] L.M. Pecora, T.L. Carroll, G.A. Johnson, and D.J. Mar, Chaos 7, 520 (1997). 
[9] G.S. Duane and J.J. Tribbia, Phys. Rev. Lett. 86, 4298 (2001). 
[10] R. Temam, Infinite Dimensional Dynamical Systems in Mechanics and Physics, Springer- 

Verlag, New- York, Applied Mathematical Sciences Series 68, (1988). 
[11] S. Weinberg, Gravitation and Cosmology: Principles and Applications of the General Theory 

of Relativity, Wiley, New York (1972). 
[12] The butterfly effect here differs essentially from bubble nucleation as in Coleman, S., Phys. 
Rev. D 15, 2929 (1977). The oscillons do not grow from a seed, but shrink, in the coordinate 
system of ([3]). Initial and final oscillon widths are both much larger than the "seed" scale. 
[13] A. Hampton and D.H. Zanette, Phys. Rev. Lett. 83, 2179 (1999). 

[14] Measure synchronization was first reported as "non-isochronic" synchronization in Ref. [3]. 
[15] J.K. Hale, J. Dynamics and Diff. Eq. 9, 1 (1997). 

[16] K. Kaneko (Ed.), Theory and Applications of Coupled Map Lattices, Wiley, New York (1993). 
[17] G.S. Duane, J.J. Tribbia, and J.B. Weiss, Nonlin. Processes in Geophys. 13, 601 (2006); G.S. 

Duane and J.J. Tribbia, Ch. 17 in Nonlinear Dynamics in the Geosciences, ed. A. Tsonis, 

Springer (2007). 



8 {May 30, 2009) 



TABLE I: Correlations between coefficients of corresponding Fourier components of the indepen- 
dent dynamical variables in the A and B subsystems of the coupled scalar field system ([5]) (a) 
and for the same system with the potential V^j^ (b). The Fourier coefficients, indexed by n, are 
partitioned into coupled modes with wavenumber k = [{n + l)/2], < n < 128 (in cycles/domain- 
length), and several ranges of uncoupled modes. Error bars for the first range of uncoupled modes 
are at 2a, where a was computed as the standard error of the mean based on a further partitioning 
into odd and even n. Uncoupled modes correlate significantly only with a nonlinear potential that 
supports oscillons. 

a) correlations with b) correlations with 

V = 02/2 - 0^/4 + 06/6 Vj^it = 02/2 + 0V4 + 06/6 
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FIG. 1: Streamfunction if) (in units of 1.48 x 10^m2s~^, averaged over layers i = 1,2) describing 
the flow at initial (a,b) and final (c,d) times, in a parallel channel model with coupling of medium 
scale modes for which \kx\ > k^o = 3 or \ky\ > kyo = 2, and < 15, for the indicated numbers n 
of time steps in a numerical integration. Parameters are as in {6^ . Synchronization occurs by the 
last time shown (c,d), despite differing initial conditions. The "blocking patterns" in the boxed 
areas coincide 

FIG. 2: Energy density p = (l/2)e~^*(02;)2 + (l/2)e^*(0t)2+e^*y(0) vs. position x for a numerical 
simulation of suggesting localized oscillons (a), and a simulation of the same equation, but 
with a different potential V{(j)) = (1/2)02 _|_ (|]^/4)^4 _|_ (^2^/5)^6^ which oscillons do not occur, 
shown for comparison (b). 
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FIG. 3: a)The local energy density p vs. x for two simulations of the oscillon system ([3]), coupled 
according to ^ and ([6|), with coupling coefficient = 2 for < 64 and = otherwise, at 
final time, {p for the second system (dashed line) is also shown reflected across the x-axis for 
ease in comparison.) The coincidence of oscillon positions is apparent.; b) Local energy density p 
vs. X at final time for two simulations of the oscillon system, plotted as in a), but with common 
initialization of all modes with k < 0.8 x 2^^^, and no subsequent coupling between the two systems. 
Oscillon positions appear uncorrelated. 



FIG. 4: a) The local energy density p vs. x for two simulations of the oscillon system ([3]), coupled 
according to ([5|) and ([6]), with coupling coefficient Cfc = 2 for k < 32 and Cfc = otherwise, at final 
time, displayed as in Fig. [31 The initial noise level was severely attenuated as compared to the 
thermal initialization used in the simulation in Fig. O the amplitude of the Fourier component 
at initial time was multiplied by (l/n)"'^^. b) Local energy density p vs. x at final time for one of 
the two simulations shown in panel a), with the part of p corresponding to < 32 shown reflected 
in the negative-?/ portion of the panel. 
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